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ABSTRACT 

We present predictions for the counts of extragalactic sources, the contributions 
to fluctuations and their angular power spectrum in each channel foreseen for the 
Planck Surveyor (formerly COBRAS/SAMBA) mission. The contribution to fluctu- 
ations due to clustering of both radio and far-IR sources is found to be generally 
small in comparison with the Poisson term; however the relative importance of the 
clustering contribution increases and may eventually become dominant if sources are 
identified and subtracted down to faint flux limits. The central Planck frequency bands 
are expected to be "clean": at high galactic latitude (\b\ > 20°), where the reduced 
galactic noise does not prevent the detection of the extragalactic signal, only a tiny 
fraction of pixels is found to be contaminated by discrete extragalactic sources. More- 
over, the "flat" angular power spectrum of fluctuations due to extragalactic sources 
substantially differs from that of primordial fluctuations; therefore, the removal of 
contaminating signals is eased even at frequencies where point sources give a sizeable 
contribution to the foreground noise. 

Key words: cosmic microwave background: anisotropies, foregrounds - radio and 
far-IR sources: counts, spatial distribution - galaxies: evolution. 



1 INTRODUCTION 

The long-sought COBE/DMR discovery (Smoot et al. 1992) 
of anisotropies of the Cosmic Microwave Background (CMB) 
has stimulated an outburst of activity in the field. Many new 
experiments have been undertaken and detections of CMB 
fluctuations have been reported on several angular scales 
(see White, Scott & Silk 1994; Smoot 1997 for reviews). 

An unavoidable fundamental limitation to these mea- 
surements is set by astrophysical foregrounds. Due to the 
large beam size, the COBE/DMR data are, to some extent, 
contaminated by the galactic emission (see, e.g., Kogut et 
al. 1996a,b) whereas they are basically unaffected by extra- 
galactic foreground sources (Banday et al. 1996; Kogut et al. 
1994). On the other hand, extragalactic foregrounds are a 
major problem for high resolution experiments reaching the 



sensitivity of AT/T ~ 10 6 , as the recently selected ESA's 
COBRAS/SAMBA - now Planck Surveyor - and NASA's 
MAP satellite missions, as well as some balloon-borne ex- 
periments, do. 

This paper presents a thorough analysis of the extra- 
galactic foreground contribution to small scale fluctuations, 
over the full wavelength range from ~ 1 cm down to ~300 
/xm covered by the Planck mission. The ten channels cur- 
rently foreseen for the experiment as well as some relevant 
details on the payload characteristics are given in Table 1. 

Estimates of temperature fluctuations due to a Poisson 
distribution of extragalactic sources have been worked out 
by Franceschini et al. (1989, 1991), Wang (1991), Blain & 
Longair (1993), Toffolatti et al. (1995), Danese et al. (1996), 
Tegmark & Efstathiou (1996), and Gawiser & Smoot (1997). 

In the frequency range of interest here, there are im- 
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portant contributions both from the high frequency tail of 
the spectrum of radio sources and from the long wavelength 
portion of the dust emission in galaxies. The present work 
improves on previous analyses for both populations. 

We adopt updated models for the evolution of galaxies 
which account for the faint number counts at 60 /jm, recently 
reassessed by Bertin, Dennefeld & Moshir (1997) as well 
as for the preliminary estimates of deep counts at 6.7 and 
15 /jm (Oliver et al. 1997) obtained from long exposures with 
the camera of the Infrared Space Observatory (ISO) and at 
170 /jm with the ISO long-wavelength photometer (Kawara 
et al. 1997), as discussed by Franceschini et al. (1997). Our 
models deal in a self-consistent way with the evolution of 
the spectral energy distribution of galaxies from UV to far- 
IR wavelengths (Mazzei, Xu & De Zotti 1992; Mazzei, De 
Zotti & Xu 1994), so that we can take advantage of sur- 
veys in a broad wavelength range to test them. Also, they 
naturally predict a relationship between the dust tempera- 
ture distribution and the intensity of star-formation activity, 
in the sense that more actively star-forming galaxies have 
warmer dust emission spectra, consistent with the observed 
correlation between the ratio /6o M m//ioo/jm and the bolo- 
metric luminosity of IRAS selected galaxies (cf. e.g. Sanders 
& Mirabel 1996). 

Moreover, we have explored various possibilities for the 
shape of the high frequency spectra of radio selected sources 
to take into account that the situation is likely rather com- 
plex: on one side, the non thermal emission spectrum of 
compact sources is expected to steepen or even to break at 
mm wavelengths; on the other hand, substantial contribu- 
tions at these wavelengths probably due to dust emission 
have been reported for several objects. 

We also give a quantitative estimate of the additional 
contributions to fluctuations due to clustering. 

The outline of the paper is as follows. In Section 2 we 
briefly review the basic formalism. In Section 3 we deal with 
source counts and their extrapolation to the Planck bands. 
Our results are presented and discussed in Section 4. Finally, 
in Section 5, we summarize our main conclusions. 

Throughout this paper we will adopt an Einstein-de 
Sitter (fi = 1) cosmology with a Hubble parameter of 
H = SOkms^Mpc- 1 . 



fluctuation level generated by randomly distributed sources 
is provided by the second moment, a, of the distribution of 
the mean number, R(x), of responses x = Sf(ip) to sources 
of flux S located at an angular distance 9 from the beam 
axis, /(i/>) being the angular power pattern of the detector, 
for which we adopt a gaussian shape with FWHM 6o- 
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where tp = (6/6o) 2 , N(S) are the differential source counts 
per steradian at a given frequency v and x c is the limit above 
which a source is considered to be individually detected. Fol- 
lowing Condon (1974), we set x c = qa; q is usually taken in 
the range 3-5. We adopt x c = 5u throughout this paper. The 
confusion standard deviation a is related to the rms bright- 
ness temperature fluctuations (AT/T rms = ((AT/T) 2 ) 1/2 ), 
at the wavelength A by: 
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is the effective beam area, dui being the solid angle dif- 
ferential element (De Zotti et al., 1996a). In terms of the 
angular distance 6 from the beam axis equation (4) reads 
Weff.i = nf d6 2 f(6) (Franceschini et al., 1989). T = 2.726 K 
(Mather et al. 1994) is the brightness temperature of the 
CMB. In terms of intensity fluctuations we have: 



/ AJ„ 



(5) 



In the Rayleigh- Jeans region AI V /I V ~ AT /T while at the 
peak of the intensity of the CMB, (A ~ 1 mm), AI V /I U ~ 
3AT/T. 



2 BASIC FORMALISM 

All the estimates presented in this paper are based on the 
assumption of "point-like" sources. As shown by Rowan- 
Robinson & Fabian (1974), this is a good approximation 
as far as the angular sizes of sources do not exceed the 
beam width. This is generally the case for the Planck mis- 
sion (FWHM> 10' at v < 140 GHz and FWHM> 4' .4 for 
the high frequency channels; see Table 1); the pixels con- 
taminated by the few bright very extended sources will be 
anyway removed. 

2.1 Sky fluctuations from randomly distributed 
sources 

The problem has been extensively discussed in the literature 
(Scheuer 1957, 1974; Condon 1974; Franceschini et al. 1989); 
we recall here only the basic points. A useful estimate of the 



2.2 Fluctuations due to clustered sources 

Clustering decreases the effective number of objects in ran- 
domly distributed cells and, consequently, enhances the cell- 
to-cell fluctuations (Peebles 1980; Barcons & Fabian 1988). 
The analysis of a complete sample of nearby (z < 0.1) 
radiogalaxies selected at 1.4 GHz (Peacock & Nicholson 
1991) has shown that, at least in a particular range of ra- 
dio power, sources are strongly clustered (correlation length 
ro — 22(ffo/50)~ 1 Mpc). Evidences of a strong angular cor- 
relation were found by Kooiman, Burns & Klypin (1995) in 
the Green Bank 4.85 GHz catalog. More recently, Loan, Wall 
& Lahav (1997) estimated the angular two-point correlation 
function of sources selected at 4.85 GHz, by combining the 
Green Bank and Parkes-MIT-NRAO surveys. For an evo- 
lution index of the correlation function (Peebles 1980) in 
the range -1.2 < e < they found 26 Mpc < (H /50)r < 
36 Mpc. 
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Table 1. Planck Surveyor — Model Payload Characteristics 



Telescope 






1.5 m Diam. Gregorian; system emissivity 1% 
Viewing direction offset 90° from spin axis 




Instrument 






LFI 










HFI 






Detector Technology 


HEMT receiver arrays 






Bolometer arrays 




Detector Temperature 






20 K 








0.1-0.15 K 






Center Frequency (GHz) 


30 


44 


70 


100 


100 


143 


217 


353 


545 


857 


Wavelength (mm) 


10 


6.8 


4.3 


3.0 


3.0 


2.1 


1.38 


0.85 


0.55 


0.35 


Bandwidth (GHz) 


6.0 


8.8 


13 


20 


37 


53 


80 


131 


202 


317 


Number of Detectors 


4 


6 


12 


34 


4 


12 


12 


6 


6 


6 


Angular Resolution 
(FWHM, arcmin) 


33 


23 


14 


10 


10.6 


7.4 


4.9 


4.5 


4.5 


4.5 


Optical Transmission 


1 


1 


1 


1 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


(AT/T) Sensitivity per res. element 
(1 year, la, 1CP 6 units) 


1.6 


2.4 


3.6 


4.3 


1.81 


2.1 


4.6 


15.0 


144. 


4630 



Significant clustering was also established for far-IR se- 
lected sources, on a scale ro ~ 8(50/-ffo) Mpc (Saunders, 
Rowan-Robinson & Lawrence, 1992) slightly smaller than 
for optically selected galaxies. 

The contribution of clustering to intensity fluctuations 
is straightforwardly obtained from the angular correlation 
function C(0*), as a function of the angular separation 
setting 6+ — (see, e.g., De Zotti et al. 1996a). If the cluster- 
ing scale is much smaller than the Hubble radius, we have: 
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where (/) is the mean background intensity and 
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Here £(r, z) = h(z)^o(r) is the two-point spatial correlation 
function, with £o(r) = (ro/r) 1 ' 8 . For the clustering evolution 
function we adopt the usual simple expression h(z) = (1 + 
z)~^ 3+e \ A(r ma x) is the value of Sz corresponding to the 
maximum scale of clustering, Si is the adopted flux limit 
and 



3eff(z) 



min[£ max ,i(S(,2)] 



dlog L L n c (L, z)K(L, z) (8) 



is the effective volume emissivity, n c (L, z) being the comov- 
ing number density of sources and K (L, z) the K-correction 
factor. The pairs of angles ip) and ($,Tp) define two di- 
rections separated by an angle #*. It is worth stressing that, 
while for differential source counts [N(S) oc S^ 13 ] with slope 
(3 < 3 the Poisson noise is dominated by the sources just be- 



low the detection threshold, Si, the main contribution due 
to clustering always comes from the faintest sources which 
do actually cluster on the given angular scale. Therefore, 
lowering Si results in an increase of the relative importance 
of the clustering term in comparison with the Poisson one. 

2.3 Angular power spectrum analysis of the 
intensity fluctuations 

We consider, as usual, the spherical harmonic expansion of 
the sky temperature fluctuations: 

fW) = E E a TYe m (^)- (9) 
e=o m=-e 

If the fluctuations are a stationary process, the angular 
power spectrum is independent of m (see Peebles 1993, p. 
517), so that 



Ct(y) 



with 



2f 



x — E <i a n 2 } = <k°i 



o° = J y(#,^(^)in, 

or, in terms of Legendre polynomials, 
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<p)\/ -^-P t (cos&) sm&dddip. (12) 



The estimate of Ct{v) is obtained by dividing the sky in 
[N x N] "equal area" cells (the cell dimension being 1.5 x 
1.5 arcmin 2 ) x [(pj-i,<fij] so that in every cell the 

fluctuation 8T/T can be considered constant. Then 



»=i j=i 
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£(») 
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where 

Pe-i (cosi?i_i) - Pt+i (costfi) + P^+1 (cOS#i-l)]. (14) 

The temperature fluctuations (ST/T) i ■ in every cell were 
calculated from simulations of the all-sky distribution of 
extragalactic point-like sources, based on our estimates of 
the source counts for each frequency channel of the Planck 
mission, rather than from the formulae derived by Tegmark 
& Efstathiou (1996). This approach allows to calculate the 
angular power spectrum under more general conditions (in 
particular allowing for the effect of clustering) since it only 
requires that the overall source distribution be statistically 
isotropic. We have checked that, assuming a Poisson distri- 
bution and using the source counts estimated by Tegmark 
& Efstathiou (1996), our method reproduces their results. 

Following Tegmark & Efstathiou (1996) we then com- 
puted the quantities 

ST t (v) = [t(t + l)C t (v)/2ir] 1/a . (15) 



3 SOURCE COUNTS AND THEIR 

EXTRAPOLATION TO MILLIMETER AND 
SUB-MM WAVELENGTHS 

3.1 Counts of radio sources 

Deep VLA surveys have allowed to extend direct determi- 
nations of radio source counts down to /ijy levels at 1.41, 
4.86 and 8.44 GHz. At these frequencies counts now cover 
about 7 orders of magnitude in flux and reach an areal den- 
sity of several sources arcmin -2 . Therefore, at these frequen- 
cies, fluctuations can be determined directly from the counts 
down to angular scales much smaller than those of interest 
here. 

As shown by Figure 1, the model by Danese et al. 
(1987) provides a good fit to the available data, at least 
for S > 100 /ijy. Particularly encouraging is the good agree- 
ment with the deep counts at 8.44 GHz (Windhorst et al. 
1993; Partridge et al. 1997), which were produced several 
years after the model, indicating that the adopted distribu- 
tion of spectral indices of sources was appropriate. 

Counts below 100 /ijy affect fluctuations due to radio 
sources on scales much smaller than those reachable by the 
Planck mission. 

The assumptions about source spectra are obviously the 
most critical ingredient for the purposed extrapolation of the 
observed counts to much higher frequencies. 

Radio loud AGNs, including "flat" -spectrum radio- 
galaxies, quasars, BL-Lacs, mostly at substantial z, are ex- 
pected to dominate the counts in the Planck low frequency 
channels for S ^ 1 mjy and, correspondingly, the contribu- 
tion of extragalactic sources to fluctuations on the angular 
scales of interest here. 

The analysis of Impey & Neugebauer (1988), giving the 
spectral indices of 162 blazars, and the quasi-simultaneous 
observations of 176 bright compact sources done by Edelson 
(1987), have shown that compact sources have "flat" spectra 
(S(u) (x v~ a , with a ~ 0, although with some scatter) at 
least up to ~ 100 GHz. This conclusion agrees with previous 



analyses (Owen, Spangler & Cotton 1980), which showed 
that the spectral indices of most strong flat-spectrum radio 
sources keep flat (a < 0.3) in the range 1 to 100 GHz. 

At still higher frequencies, a steepening or even a spec- 
tral break of the synchrotron emission is expected. On the 
other hand, the presence of excess mm emission over ex- 
trapolations from cm wavelengths has been established by 
Knapp & Patten (1991) for a sample of nearby radio galax- 
ies. Observations of large mm fluxes attributed to dust emis- 
sions have been reported for several distant radio galaxies 
(see Mazzei & De Zotti, 1996 and references therein). The 
inferred dust masses are 1-2 orders of magnitude higher than 
found for nearby radio galaxies. The two components (syn- 
chrotron and dust emission) may well have different evolu- 
tion properties. 

In view of the uncertainties on the spectra mentioned 
above, our estimates of counts and fluctuations due to radio 
selected sources have been calculated as it follows. We adopt 
the Danese et al. (1987) simple luminosity evolution model 
and three different choices for the average spectral index of 
"flat-spectrum" compact sources: a) a = —0.3, b) a = 0.0, 
and c) a = 0.3 at 20 £ v ;$ 200 GHz, with a steepening to 

a = 0.7 at higher frequencies; below 20 GHz we have set 
a = 0. As for "steep"-spectrum sources (elliptical, SO and 
starburst galaxies), whose contribution to source counts is 
actually minor in the whole frequency range of interest here, 
the radio power - spectral index relation determined by Pea- 
cock and Gull (1981) has been adopted. The maximum and 
minimum number of expected radio sources in each Planck 
channel, given in Table 2, correspond to the two extreme 
values of a for compact sources. 

Holdaway et al. (1994) carried out sensitive 90 GHz ob- 
servations of a sample selected at 5 GHz and observed at 8.4 
GHz. They derived a distribution of spectral indices between 
8.4 GHz and 90 GHz which allowed them to extrapolate to 
90 GHz the 5 GHz source counts. They estimate that, over 
the entire sky, there are 178 sources with Sgo GHz > 1 Jy, al- 
most a factor of 2 less than predicted by our model. However, 
as discussed by these authors, their estimate is somewhat 
below the number of known sources with Sgo GHz > 1 Jy, so 
that it should be viewed as a lower limit. 

On the other hand, our estimate may miss a population 
of sources with strongly inverted radio spectra. However, the 
analysis by Condon et al. (1995) of a large sample of extra- 
galactic sources detected by IRAS at A = 60 pm and identi- 
fied with VLA radio source catalogues at 4.85 GHz, found no 
evidence of a significant population of sources with spectra 
rising steeply from centimeter to millimeter wavelengths. 



3.2 Counts of far— IR sources 

Both in the case of normal and of many active galaxies (as 
far as we can tell, based on the very limited information cur- 
rently available) , at wavelengths shorter than a few mm (in 
the rest frame) , dust emission rapidly overwhelms the radio 
emission. Due to the very steep increase with frequency of 
the dust emission spectrum at mm and sub-mm wavelengths 
(q ~ —3.5) the wavelength at which dust emission takes over 
does not change much between radio quiet and radio loud 
sources, in spite of the fact that the ratio of radio to far-IR 
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emission for the former class is orders of magnitude lower 
than for the latter. 

As in the case of radio sources, estimates of the counts 
of far-IR sources in the frequency region covered by the 
Planck instruments are difficult because of the wide gap with 
the nearest wavelength (60 (im) where the most extensive 
(yet relatively shallow) surveys exist. There is a considerable 
spread in the distribution of dust temperatures, so that the 
observed 1.3 mm/60 /im flux ratios of galaxies span about 
a factor of 10 (Chini et al. 1995; Franceschini & Andreani 
1995; see also De Zotti et al. 1996b). A tentative estimate of 
the luminosity function of galaxies at mm wavelengths based 
on a 1.25 mm/60 /im bivariatc luminosity distribution has 
been presented by Franceschini, Andreani & Danese (1997). 

Furthermore, the observational constraints on evolution 
properties of far-IR sources are very poor. They come essen- 
tially from IRAS 60 /jm counts, which cover a limited range 
in flux and are rather uncertain at the faint end (Hacking & 
Houck 1987; Gregorich et al. 1995; Ashby et al. 1996; Bertin, 
Dennefeld & Moshir 1997). 

From a theoretical point of view, there is a great deal 
of uncertainty on the physical processes governing galaxy 
formation and evolution. Some models assume that the co- 
moving density of galaxies remained essentially constant af- 
ter their formation, while they evolved in luminosity due to 
the ageing of stellar populations and the birth of new gener- 
ations of stars (pure luminosity evolution; see, i.e., Frances- 
chini et al., 1994 for a most complete discussion on the 
subject). On the other hand, according to the hierarchical 
galaxy formation paradigm, big galaxies are formed by coa- 
lescence of large numbers of smaller objects (see, e.g. White 
1996 and references therein). 

Furthermore, evolution depends on an impressive num- 
ber of unknown or poorly known parameters: merging rate, 
star formation rate, initial mass function, galactic winds, 
infall, interactions, dust properties, etc. 

Although the evolutionary history is highly uncertain, 
strong evolution is expected in the far-IR/mm region partic- 
ularly for early type galaxies since during their early phases 
they must have possessed a substantial metal enriched inter- 
stellar medium. This expectation is supported by evidences 
of large amounts of dust at high redshifts (cf. in partic- 
ular the case of high- 2: radiogalaxies: Mazzei & De Zotti, 
1996) and by the intensity of the isotropic sub-mm com- 
ponent reported by Puget et al. (1996) from an analysis of 
COBE/FIRAS data (see Burigana et al. 1997). 

Moreover, the K-correction due to the steep increase 
with frequency of the dust emission spectrum below the 
peak, generally located at A ~ 60-150 fim, strongly ampli- 
fies the evolutionary effects, that may thus be appreciable in 
the relatively shallow Planck surveys, at least in the highest 
frequency bands. 

In view of the above, we adopt here two different de- 
scriptions of the cosmological evolutions of far-IR sources. 

i) Our reference model is an updated version of model 
C of Franceschini et al. (1994; opaque model), which pro- 
vides a good fit (see Fig. 2) to the 60/im IRAS counts, as 
recently reassessed by Bertin et al. (1997) and to the far-IR 
extragalactic background spectrum estimated by Puget et 
al. (1996), as shown by Burigana et al. (1997). The model 
also predicts counts at 170 /im and at 850 fim which are in 
agreement, respectively, with the preliminary estimates by 



Kawara et al. (1997) using the ISO long-wavelength pho- 
tometer and by Smail, Ivison & Blain (1997) using the new 
bolometer array, SCUBA, on the JCMT (see Fig. 3). 

Moreover, we find agreement with the recent estimate of 
the "source" confusion noise by Bertin et al. (1997; their Sec- 
tion 4.1), from the analysis of the IRAS Very Faint Source 
Sample (VFSS). Their estimated value is CT con f ~ 20 mjy, 
about 2/3 of the total noise of 30.2 ± 1.2 mjy (instrumen- 
tal noise plus source confusion added up in quadrature). By 
comparison, our reference model gives a confusion standard 
deviation of cr con f ~ 18 — 20 mjy, with a gaussian beam pat- 
tern of FWHM~ 3 arcmin and integrating up to a source 
detection limit x c ~ 100 mjy. 

The updated opaque model adds a density evolution up 
to z = 2 of late type galaxies (spirals, irregulars, starburst) 
only, proportional exp(2r(z)), with t(z) = 1 — t(z)/to, t(z) 
being the age of the universe at redshift z and to its present 
age. The estimated AT/T rms levels as well as the spatial 
power spectra of fluctuations due to far-IR extragalactic 
sources have all been computed exploiting this model (see 
Section 4). 

On the other hand, the signal detected by Puget et al. 
(1996) may be contaminated by cold dust emission from an 
extended Galactic halo (Fixsen et al. 1996). In that case 
our reference model would provide an upper limit to the 
evolution of extragalactic sources in the far-IR. 

it) A lower evolution rate in the far-IR is implied by 
the moderate extinction model by Franceschini et al. (1994). 
This second model predicts source counts which are still 
compatible with the IRAS 60 fim counts and with the quoted 
estimated counts at 170 /im by Kawara et al. (1997). On 
the other hand, they slightly underestimate the number of 
sources quoted by Smail, Ivison & Blain (1997). 

The estimated number counts, N(> Shu^fir, of far-IR 
sources in the HFI frequency channels are a factor ~ 1.2 — 1.8 
lower than the counts plotted in Fig. 4 (our reference model) 
if Siim — 1 Jy (the lower value, 1.2, applies to the number 
counts ratio at 857 GHz whereas the value 1.8 is found for 
the same ratio at 143 GHz). If we adopt the fainter flux limit 
of Sum ^ 10~ 2 Jy then the model predicts counts ~ 1.4 — 2.2 
times lower than the updated model C. The increasing dif- 
ference between the estimated number counts, if one moves 
shortward in frequency from 857 GHz down to 143 GHz, 
is determined by the steep rise with increasing frequency of 
the dust emission spectrum in the sub-mm domain which en- 
tails more important contributions to the counts from high 
redshift galaxies. 

The AT/Trms levels are only ~20-40% lower than those 
predicted by the reference model at the angular resolution 
limit of the HFI instrument. At larger angular scales the 
difference is negligible since the two models predict a very 
similar number of bright sources (see before). The maximum 
and minimum number of expected individually detected far- 
IR sources in the Planck surveys at \b\ > 20° given in Table 
2 correspond to the predictions of the updated model C and 
of the moderate extinction model. 

The contribution of dust-rich galaxies to the radio 
counts (see Table 2) has been estimated assuming a linear 
relationship between the far-IR flux at 60 /u,m (dominated 
by the starburst component) and the radio centimetric flux. 
The ratio of the 60 /j,m to the radio flux at 1.4 GHz was as- 
sumed to be S60/jm/Si.4GHz — 140 (see Helou et al., 1985). 
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The contributions of AGNs to the counts in the high 
frequency channels (where, however, star-forming galaxies 
should be the dominant population) are even more uncer- 
tain, since they depend on the evolution of both the non- 
thermal component and the dust emission, which in turn 
depend on several unknown or poorly known factors, such 
as the evolution of the nuclear energy source, the effect of 
possible circumnuclear starbursts, the abundance, proper- 
ties and distribution of dust, and so on. The IRAS survey 
data do not help much, since the detection rate of quasars 
was extremely low. 

Taking only into account the evolution of the non- 
thermal component, assumed to parallel that observed in 
X-rays, and adopting the spectral energy distributions used 
by Granato et al. (1997) for type 1 and type 2 AGNs, we 
expect a detection rate of radio quiet AGNs increasing with 
increasing frequency: only a few of them are expected at 1.4 
mm, but their number should increase up to a few hundreds 
at the highest frequencies. In any case, the number of ra- 
dio quiet AGNs is negligible compared with the number of 
far-IR galaxies. 

3.3 Total counts 

Our estimates (see Table 2 and Figure 4) indicate that 
Planck counts are dominated by radiosources for frequen- 
cies up to about 200 GHz. Note that the large differences 
in the estimated numbers of radio selected sources is deter- 
mined by the variation of the adopted values for the average 
spectral index of compact sources (see Section 3.1). It is in- 
teresting that significant numbers of radio galaxies should 
be detectable in all frequency channels, so that the Planck 
mission will allow to explore their poorly known high fre- 
quency spectra. 

Also, the number of galaxies that should be detected 
in the high frequency channels is large enough to allow sta- 
tistical investigations of their properties in this particularly 
interesting range and the definition of reliable local luminos- 
ity functions. 



4 RESULTS AND DISCUSSION 

Our estimates of the rms temperature fluctuations 
((AT/T) rms ee ((AT/T) 2 ) 1/2 ) as a function of the angu- 
lar scale, and of the angular power spectra [5Ti, eq. (15)] of 
fluctuations due to extragalactic point sources are shown in 
Figures 5-6 and 7-8 respectively. 

4.1 Poisson fluctuations 

In Figures 5 and 6 the thick solid line corresponds to the case 
of "flat" -spectrum radio sources with average spectral index 
a — above 20 GHz and far-IR sources evolving accord- 
ing to the updated model C by Franceschini et al. (1994). 
Sources brighter than 1 Jy are assumed to be individually 
identified and removed (a rather conservative assumption: 
cf. Table 2). The effect of subtracting fainter sources, identi- 
fied by means of independent, deeper surveys is also shown. 

Although we expect that the individually detected 
sources in the high galactic latitude Planck surveys at fre- 
quencies up to about 200 GHz (see Table 2) will be mostly 



radio sources, a significant contribution to fluctuations may 
come from far-IR sources well below the detection limit. 

In fact, the steepness of the counts of far-IR sources 
in the flux density interval where evolution sets in (and 
its effect is boosted by the steep K-correction) implies that 
their dominant contribution to fluctuations may not come, 
as usual, from flux densities corresponding to one source per 
beam (where radiosources probably dominate) but, in some 
cases, from the fainter fluxes where the slope /3 of differential 
counts becomes > 3. 

On an angular scale 9o — 10' far-IR sources do not con- 
tribute more than a few % of the confusion fluctuations for 
v £ 100 GHz; at 150-200 GHz the contributions of radio and 
far-IR sources are comparable, while at higher frequencies, 
far-IR sources dominate. 

The amplitude of confusion fluctuations, after subtrac- 
tion of only those sources which are individually detected by 
the Planck instruments themselves, decreases from AT/T ~ 
5 x 10~ 6 at 30 GHz to ~ 10~ 6 at ~ 100-150 GHz, well below 
the expected rms primordial fluctuations. 

It may be noted that up to ~ 100 GHz the rms fluctu- 
ations due to discrete sources, as a function of the angular 
scale, show a broad maximum at #o ~ 10', while they in- 
crease with decreasing 9o at higher frequencies, due to the 
contribution of evolving far-IR sources which becomes more 
and more important with increasing angular resolution of 
the survey. 

As previously pointed out, the estimates of fluctuations 
due to far-IR sources are particularly uncertain since the 
available information on both their spectra and their evolu- 
tionary properties is poor. However, the deep IRAS 60 /jm 
counts on one side and the upper limits on the far-IR back- 
ground intensity on the other, strongly constrain the am- 
plitude of fluctuations at the angular resolution of the high 
frequency channels of the Planck mission. Phenomenologi- 
cal pure luminosity evolution models without changes in the 
source spectra with cosmic time (Franceschini et al., 1988; 
1991) consistent with the high 60 fim counts by Gregorich 
et al. (1995), which may be overestimated at the faint end 
because of source confusion, and barely consistent with the 
conservative upper limit on the extragalactic far-IR back- 
ground derived by Shafer et al. (1997), entail fluctuations 
very close to those implied by the updated model C on scales 
above the angular resolution of the Planck mission; only be- 
low ~1.5 arcmin the fluctuation level increases by a factor 
of 1.5-2. 

Note that the effect of source variability is already in- 
cluded, in a statistical sense, in the above estimates. In fact, 
variability affects source counts in a manner similar to the 
Eddington effect: the observed counts are made sistemati- 
cally higher at brighter fluxes, since sources are preferen- 
tially detected in their brighter phases. Since fluctuations 
are estimated from observed counts, the effect is automati- 
cally taken into account. On the other hand, variability seri- 
ously hinders subtractions of sources based on observations 
at different frequencies or non-simultaneous. 



4.2 The effect of clustering 

In the case of radio sources, based on the results by Loan 
et al. (1997), we have adopted an angular two-point correla- 
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Table 2. Instrumental and confusion noise estimates for the Planck mission and expected numbers of individually 
detectable sources at \b\ > 20°. 



^cff Kb beam cr noiso <r Gal cr conf - °" CMB S lim N(> 5i im ) rad io N(> Sii m )FlR 

(GHz) (mm) (arcmin) (mjy) (mjy) (mjy) (mjy) (mjy) (8 sr) (8 sr) 



30 


10.0 


33 


13 


100 


37 


78 


650 


300-800 


0-2 


44 


6.8 


23 


19 


45 


17 


80 


480 


450-1600 


3-10 


70 


4.3 


14 


25 


15 


8 


68 


330 


600-2500 


6-18 


100 


3.0 


10. 


27 


7 


6 


64 


350 


500-2800 


4-12 


100 


3.0 


10.6 


13 


7 


6 


64 


330 


500-2800 


4-12 


143 


2.1 


7.4 


12 


6 


4 


56 


290 


330-3600 


5-15 


217 


1.38 


4.9 


14 


5 


4 


31 


180 


350-4000 


120-200 


353 


0.85 


4.5 


2d 


18 


18 


16 


200 


250-2000 


2000-3500 


545 


0.55 


4.5 


44 


62 


45 


3 


450 


40-250 


6000-10000 


857 


0.35 


4.5 


36 


120 


70 




700 


10-50 


17000-25000 



f 1 ) Flux density fluctuations of the cosmic microwave background corresponding to AT/T = 10 5 

( 2 ) Adopted detection limit for discrete sources, equal to 5 X (o'^ lolsc + CT conf + ""Galaxy "cmb) 1 ' 2, ^ rough estimate of anisotropies, 
CT Galaxyi due to the galactic emission at high galactic latitudes > 20°), has been obtained following Danese et al. (1996) and Bcrsanclli 
et al. (1996; Fig. 2.3, 2.4 and Appendix A.l) and taking into account the angular dependence of galactic dust fluctuations determined 
by Gautier et al. (1992). If we neglect this correction factor, the galactic dust fluctuations at 217, 353, 545 and 857 GHz turn out to be 
larger by a factor ~ 2.2. 



tion function w(8) = 0.01 for 9 = 1°. For the double power 
law representation of £(r, z), described in Section 2.2, this 
implies a clustering radius ro = 36(50/.ffo) Mpc if e = 
(stable clustering); if e = —1.2 (clustering constant in co- 
moving coordinates), ro = 26(50/ Ho) Mpc. The two models 
yield differences on the estimated contribution of clustering 
to fluctuations of 10-15%. The results shown in Figs. 5 and 
6 refer toe = —1.2. 

For far-IR sources we have adopted r = 10 Mpc (Ho = 
50) and e — 0. This value of the clustering radius is slightly 
larger than that estimated by Saunders et al. (1992) for far- 
IR selected sources, on account of the fact that, according 
to our reference model, a substantial contribution to fluc- 
tuations comes from early phases of the evolution of early 
type galaxies, which are more clustered than the late type 
galaxies studied by Saunders et al. (1992). 

The adopted evolutionary and spectral properties of the 
sources are the same as for the calculations of Poisson fluc- 
tuations. In all cases the latter dominate except in the case 
that sources can be identified and subtracted down to fluxes 
below the detection limit of the survey, thus greatly decreas- 
ing the Poisson fluctuations, while the contribution due to 
clustering is only weakly affected. 

As shown by Fig. 6, the autocorrelation of intensity 
fluctuations essentially reflects the clustering of radiosources 
below 200 GHz and of far-IR sources at higher frequencies; 
around 200 GHz the contributions of the two populations 
are similar. 



4.3 Angular power spectra of fluctuations 

Following Tegmark & Efstathiou (1996) we have computed, 
and plotted in Figures 7 and 8, the quantity 5Ti(v) = 
[£(£ + l)Ce{v) /2-k] 1 ^ 2 which is roughly the average rms tem- 



perature fluctuation in the multipole range £ < £ < £\ with 
ln(V4>) = L 

The plotted power spectra are the mean of 50 simu- 
lations of the all sky distribution of sources. Fluctuations 
are dominated by radio sources at the lowest frequencies 
(30, 44, 65 and 100 GHz) and by the far-IR sources at the 
highest (353, 545 and 857 GHz); both classes of sources 
are important at the intermediate frequencies (143, and 217 
GHz). For far-IR sources the updated model C has been 
adopted. As mentioned above, this model may overestimate 
the fluctuation level. A spectral index a — was adopted 
for "flat-spectrum" radiosources. The effect of clustering was 
neglected. 

It is assumed that sources brighter than 1, 0.1 and 0.01 
Jy are identified and removed from the maps (see caption). 
For the 217 GHz band we also show the effect of removing 
only the few brightest sources (S > 10 Jy). As shown by 
Table 2, we expect that it will be possible to directly detect, 
and subtract out, sources down to flux densities of a few to 
several hundred mjy. 

It may be noted that lowering the flux limit for source 
subtraction decreases the amplitude of fluctuations more ef- 
fectively at low than at high frequencies. This is because 
of the steepness of the counts of far-IR sources at sub-mm 
wavelengths, implying substantial contributions to fluctua- 
tions from very faint sources. 

Our results for a flux density cutoff of 100 mjy are 
substantially below those obtained by Tegmark & Efstathiou 
(1996) adopting the same flux limit (see Fig. 7). This is due 
to their assumption of a spectral index a = for all sources 
selected at 1.5 GHz. This leads to a substantial overestimate 
of the expected counts, since most 1.5 GHz sources in the 
relevant flux density range are known to have a steep (a ~ 
0.7) spectrum. 
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5 CONCLUSIONS 

Although, as stressed above, the present estimates are rather 
uncertain, particularly in the high frequency bands, it can be 
safely concluded that at least the central frequency channels 
will allow a clean view of primordial anisotropics up to the 
maximum ^-values accessible to the Planck mission. In fact, 
even under the most extreme assumptions compatible with 
the available information on high frequency spectra of radio 
sources and on the evolution properties of far-IR sources, 
the amplitude of fluctuations due to discrete sources in the 
100-200 GHz range are well below the expected amplitude 
of primordial anisotropies. 

The availability of multifrequency data allows an ef- 
ficient identification of pixels contaminated by discrete 
sources. Table 2 shows that many sources not directly de- 
tectable in the central frequency channels can be identified 
at higher or lower frequencies. The corresponding pixels can 
be simply removed; thanks to the large area surveyed, the 
number of remaining clean pixels will be nevertheless very 
high. But a reliable subtraction of the contaminating flux 
may also be possible thanks to the fact that a significant 
number of sources of the various classes should be detectable 
in most or all frequency channels, allowing the definition of 
well defined template spectra. 

Moreover, removal of contaminating signals is eased 
by the substantial difference between their power spectrum 
and that of primordial fluctuations (Tegmark & Efstathiou 
1996). 

On the other hand, the Planck mission will bridge the 
gap between radio surveys, carried out at v < 8.4 GHz (A > 
3.6 cm) and far-IR surveys (IRAS and ISO) at v > 1500- 
3000 GHz (A < 100-200 /im), providing the first exploration 
of the whole sky in a spectral region where many interest- 
ing astrophysical phenomena are most easily investigated. A 
pot-pourri of issues for which these data will be extremely 
relevant include: bremsstrahlung emission as a tracer of evo- 
lution of stellar populations; high-frequency behaviour of 
the spectra of compact radio sources and implications for 
their physical properties; definition of unbiased samples of 
blazars; cold dust in galaxies and hints on its evolution; 
physical and evolutionary connections between nuclear ac- 
tivity and processes governing the abundance and the prop- 
erties of the interstellar material; the relationships between 
different AGN classes and tests for unified models; energy 
source(s) of the huge far-IR emission from type 2 Seyferts 
and from some QSOs and radiogalaxies. 
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FIGURE CAPTIONS 

Figure 1. Comparison between predicted and observed 
differential source counts at 1.4, 5 and 8.44 GHz normalized 
to 150S~ 2 ' 5 sr -1 Jy -1 (see top panel). The contributions of 
the most relevant classes of radio sources according to the 
model of Danese et al. (1987) are shown. For references for 
the data points see Danese et al. (1987); additional data are 
from Donnelly et al. (1987), Fomalont et al. (1991), Wind- 
horst et al. (1993). 

Figure 2. Differential 60 /im galaxy counts normal- 
ized to 600S^ 2 - 5 sr _1 Jy -1 . Data are from Hacking & Houck 
(1987), Rowan-Robinson et al. (1991), Gregorich et al. 
(1995), and Bertin et al. (1997). The solid line shows the pre- 
dictions of model C by Franceschini et al. (1994), updated 
as described in the text. Also shown are the contributions 
to the counts from late type galaxies (spiral+irregular, short 
dashed line; starburst, dotted line), from early type galax- 
ies (E and SO), assumed to go through a strongly dust ab- 
sorbed phase during their early evolution (dot-dashed line), 
and from AGNs (long/short dashes). 

Figure 3. Integral counts at 170 fim and at 850 fim. The 
meaning of the lines is the same as in Figure 2. The data 
points are the preliminary estimates by Kawara et al. (1997) 
at 170 fim and by Smail et al. (1997) at 850 fim. As stressed 
by the latter authors, their result might be, conservatively, 
taken as an upper limit; the plotted error bars for this point 
are 2a. 

Figure 4. Predicted integral counts in the Planck 
bands. The solid line shows the total counts. The dashed 
and dotted lines show, respectively, the contributions of ra- 
dio and far-IR selected sources, as predicted by the models 
by Danese et al. (1987) and the Franceschini et al. (1994), 
updated as described in the text. We have adopted an aver- 
age spectral index a = for "fiat" -spectrum radio sources 
and the "opaque" model for the early evolution of early-type 
galaxies in the far-IR. 

Figure 5. Temperature fluctuations due to discrete 
sources at i/=100 GHz as a function of the angular scale. 
The solid, short-dashed, long-dashed, and dot-dashed lines 
covering the full range of angular scales correspond to the 
estimated Poisson contributions for different choices of the 
limiting flux above which sources are individually detected 
and subtracted out: 1 Jy, 0.1 Jy, 10 mjy, and 1 mjy, respec- 
tively. The lines plotted for scales larger than the Planck res- 
olution at this frequency (10') show the contributions due to 
clustering of radiosources (clustering of far-IR sources has a 
negligible effect at this frequency); the upper curve assumes 
that sources are subtracted for S > 1 Jy, the lower one, for 
S > 0.1 Jy. 

Figure 6. Temperature fluctuations due to discrete ex- 
tragalactic sources at the frequencies indicated in each panel, 
as a function of the angular scale. As for the Poisson con- 
tribution, the lines have the same meaning as in Figure 5. 
Again, the contributions due to clustering are shown only for 
scales larger than the angular resolution of Planck instru- 
ments at each frequency: the dotted and dot-dashed lines 
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correspond to radiosources and to a flux limit of 1 and 0.1 
Jy for source subtraction, respectively; the solid line corre- 
sponds to far-IR sources and a to S; = l Jy for source sub- 
traction (lowering S; down to 0.1 Jy provides a negligible 
variation of the estimated contribution). At frequencies up 
to ~ 100 GHz the effect of clustering of far-IR sources is neg- 
ligible and only the contribution of radiosources is shown; 
the opposite is true at v > 353 GHz. The contributions of 
both populations are given at 143 and 217 GHz; at the latter 
frequency they are very close to each other. 

Figure 7. Angular power spectra of fluctuations due 
to a Poisson distribution of discrete sources at v = 100 GHz 
(see text for details). The solid, dot-dashed, and long/short 
dashed lines correspond to a flux limit for source removal 
of 1, 0.1, and 0.01 Jy, respectively. The dotted lines show 
the "unsmoothed" noise fields (defined as in Tegmark & Ef- 
stathiou 1996) foreseen for the Planck Surveyor detectors. 
The roughly horizontal thin dot-dashed line shows the de- 
pendence on multipole of temperature fluctuations predicted 
by the standard CDM model (scale-invariant scalar fluctu- 
ations in a Q, — 1 universe with Ho = 50kms _1 Mpc -1 , 
and baryon density Qb = 0.05). The thin long-dashed line, 
labelled T&E96, shows the power spectrum estimated by 
Tegmark & Efstathiou (1996). 

Figure 8. Same as in Figure 7, for the frequencies 
indicated in each panel. The dashed line in the 217 GHz 
panel refers to a removal of only the few brightest sources 
(S > 10 Jy). The short-dashed line in the upper left cor- 
ner of the 30 GHz panel shows, for comparison, the "un- 
smoothed" noise field for COBE/DMR, calculated adopting 
the average pixel noise given by Bennett et al. (1994). 
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